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Abstract 

We consider a general spin- 1/2 Ising model with multisite interaction on the 
Husimi lattice with the coordination number q and derive an analytical expression of 
correlation functions for stable fixed points of the corresponding recurrence relation. 
We show that for q = 2 our model transforms to the two-state vertex model on the 
Bethe lattice with q = 3 and for the case q = 3, with only nearest neighbour 
interactions, we transform our model to the corresponding model on the Bethe 
lattice with q = 3 , using the Yang-Baxter equations. 
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1 Introduction 



Correlation functions play an important role in the understanding of the critical behavior 
of the statistical mechanical models. There were received some exact results for corre- 
lation functions and critical exponents for the Euclidean plane models using conformal 
invariance |1|, 0. Usually Bethe-like (Husimi) lattices considered as embedded in the 
infinite-dimensional Euclidean space. However Bethe-like lattices can be embedded in a 
two-dimensional space of constant negative curvature (the hyperbolic plane) with fixed 
bond angles and lengths. These hierarchical structures like the Bethe and Husimi lattices 



can be described as a conformal (regular) tiling of the hyperbolic plane [10, 11]. It is 
believed that several among its interesting thermal properties could persist for regular 
lattices, for which the exact calculation is so far intractable. 

The exact correlation function for Ising model on the Bethe lattice was obtained by 
Mukamel || and Falk [0] and for Potts model on the Cayley tree in zero magnetic field 
by Wang and Wu ||. The two- vertex model on the Bethe lattice with the coordination 
number q = 3 was explored [Q, a closed formula for correlation functions in stable fixed 
points of the corresponding recurrence relation was derived and two types of first-order 
phase transitions which were distinguished according to the behavior of correlation func- 
tions in the high-temperature phase were obtained || . Recently Izmailian and Hu derived 
the exact correlation functions for the Ising, BEG and the most general spin -S models 
on the Bethe lattice j|, |5|. Applying combinatorial approach and graph expansion Yang 
and Xu calculated multi-spin correlation function of the Ising model on Bethe-type lattice 
M, they results indicates no existence of long range correlation at any finite temperature 
which implies that no finite temperature phase transition occurs. Recently, the multi- 
site interaction Ising model on the Husimi lattice was investigated. First, it was shown 



that this approach yields good approximation for the ferromagnetic phase diagrams [|14 



and closely matches the exact results obtained on a Kagome lattice [0. The reason for 
studying the Ising model with multisite interaction is that it plays an important role in 
the investigations of real physical systems such as binary alloys |17|, classical fluids [18] , 
solid 3 He [TIJ, liquid bilayers PPJ , and rare gases |2~T| . 

In this paper we derive a closed formula for correlation functions of the multisite 
interaction Ising model on the Husimi lattice for arbitrary coordination number q and 
show that our model for the case q = 2 coincides with the two-vertex model and for 
only nn (nearest neighbour) interaction Ising model it transforms into the corresponding 
model on the Bethe lattice with q = 3. 

The plan of the paper is as follows. In section 2, we define the model, derive the 
relevant recurrence relation and magnetization deep within the tree. In section 3, we 
derive a closed formula for correlation functions g{n) = (SiS n ) — (Si) (S n ) , where the 
symbol (. . .) signifies the mean value and Si and S n are spins on 0th and {n — l)th 
shells respectively. Section 4 is devoted to the star-triangle relation between the partition 
functions of Ising models on the Husimi and Bethe lattices. 



2 Recurrence relation and magnetization 



The pure Husimi lattice [IT2L Of, shown in Fig.l, is characterized by the coordination 



number q, the number of triangles which goes out from each site, the Oth-generation is a 
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single central triangle. 




Fig.l The Husimi lattice with q = 3 . 

The general Hamiltonian with pair, three-site interactions in an external magnetic 
field is 



K1 (ij) A 



(1) 



where S% takes values ±1; the first sum goes over all nearest neighbour pairs of the lattice 
sites, the second sum goes over all triangular faces of the Husimi lattice and the third 
one goes over all sites. The advantage of the Husimi lattice introduced here is that for 
the models formulated on it, the exact recurrence relation can be derived. From Fig.l 
it is apparent that if the graph is cut at the central triangle, then it splits up into three 
identical disconnected pieces and each of them contains 7 sub-lattices, where we defined 
7 = 9—1. Then the partition function 



Z(N) 



{S} 



H 

kf 



may be written 
Z(N) 



£ exp 

{So} 



h 



J 2 (S S -\- S Sq -\- S Sq) + J^SqSqSq + ~i^o ^( 



x (2) 



xB\N,Sti&(N,S*)&(N,S$) 



where Sq are spins of the central triangle and B(N, S) is the partition function of one 
sub-lattice containing iV shells and with the root-spin fixed in the state S . Here we use 
the same notations as in Kolesik paper 0. Because of our Husimi lattice is constructed 
from triangles we can introduce statistical weights depending on three spins located on 
triangle's vertexes. Introducing the notations 



w(S 1 ,S 2 ,S 3 ) = exp 



MSiS 2 + S 1 S 3 + S 2 S 3 ) + J 3 S 1 S 2 S 3 + -(S x + S 2 + S 3 ) 

Q 



(3) 
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and 
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a-2 



a 3 
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= exp 
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+ 3- 



-J 2 + J 3 - - 

-J 2 - J 3 



we can present the partition function as follows 

Z(N)= £ ^(5 1 ,5 2 ,5 3 )^(iV,5 1 )^(iV,5 2 )^(iV,5 3 ; 

ol q2 q3 



(4) 
(5) 
(6) 
(7) 

(8) 



Using the "cutting apart" procedure described above, one can derive the recurrence rela- 
tion for B(N, S) 



B(N, S)= J2 W ( S > Si, S 2 )Bi(N - 1, SJB^N - 1, S 2 ). 

Si,S 2 

The magnetization of the central site we can express as follows 



(9) 



m 



So 1 ) = (Z(N))- 1 £ S l w (Si Si Si) Br(N, S l )B^(N, S 2 )B^(N, S t 

cl c2 q3 
°0' D 0'°0 



Let the following variables be introduced 

q N {S) - 



and xn 



B(N, -) 



B(N,+) " ~" B(N,+Y 
Then using (^) one can write the recurrence relation for xn 

a 2 + 2aiX 1 + a x 21 



(10) 



x N = f (x N -i) , f(x) 



a 3 + 2a 2 x 7 + aix 27 



The xn has no direct physical meaning, but through it one can express the magnetization 
of the central site 



a 3 + a 2 Xjy — aix 2 ^ — aoX N ' 
m " a 3 + 3a 2 ^ + 3ai 4 7 + ao4 7 ( j 

and other thermodynamic parameters, since we can say that the xn determine the states 
of the system. 

In this paper we restrict our treatment to ferromagnetic case ( J 2 > and J 3 > 0) 
and antiferromagnetic one (J 3 < or/and J 2 < 0) at high temperatures when our one 
dimensional mapping ( [TT]) has stable fixed points only. For the antiferromagnetic case at 
low temperatures T < T* the single fixed point Xq become unstable when 

~Q~^ f( X ) \x=X < —1 



Jh 
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and a so-called period doubling bifurcation occurs: the recursive sequence {xjv} converges 
now not to the single fixed point but to the stable two-cycle {cci,^}- This phase should 
be explained as an arising of a two-sublattice phase such that x\ and x 2 determine the 
states on each sublattice JT2], [13], In terms of one dimensional maps the condition of 



the second order phase transition point of ferromagnetic models is 

f(x) \ x=x = 1 

d x 



and the second order phase transition point of antiferromagnetic models 15 



is 



dx f(x)\ x=x - 1. 

We restrict our treatment to stable fixed points of the mapping (|Tl|) because the procedure 
applied later usable for fixed points only (see Sect. 3). The fixed point x of the recurrence 
relation (|Tl|) is a solution of the equation 

x = f (x) or a\x 2l+l — a^x 21 + 2a 2 x 1+1 — 2a\x 1 + a^x — a 2 = 0. (13) 

Thus, for the stable fixed point x one obtains for the magnetization 

r 1 — ^at^n+i I — x 9 , . 

m = hm 7z = . (14) 

JV^ool + X J N X N+1 1 + Xl 

3 Correlation Functions 

One can define correlation functions on the Husimi lattice as follows 

g(ri) = (SiS n ) - (Si) (S n ) 

where spin Si belongs to the central triangle (0th shell) and S n located somewhere in the 
{n — l)th shell. Using denotation ([U]) and recurrence relation @ we can express (SiS n ) 
as follows 



(SiS n ) = lim 



E SiqJf(Si)w(Si,yi, zi)qh(yi)qk 1 (zi)w(zi, y 2 , z 2 )q 1 N _ 1 (y 



2, 



N ^°° EqN(Si)w(Si, yi, z^qjfiy^qjf 1 (z l )w(z 1 , y 2 , ^2)^-1(2/2 

QN-n+2( Z n-l) w ( z n~li Un, S n )qJ{ (2/ n )<?Ar -n+li^n) S n 



(15) 



' ' ' QN-n+2( Z n-l) w { z n~l, yni ^n)?N-n+l (^n)?AT-ra+l (^n) 

Let us define following matrices with elements 

M(zi,z i+ i) = ^q]^} i+ i(z i )w(z il y i+u z i+ i)q r N _ i (y i+ i), (16) 



Vi+l 



A'(S n ,Si) = S n Siq N (Si)qJ lf _ n+ i(S n ), (17) 
A(S n ,Si) = qN(Si)q]f„ n+ i(S n ). (18) 

Note that for stable fixed points of the mapping ([TTJ) the difference between qN(S) and 
QN-i(S) disappears for TV — > 00 and for all finite values % . Hence, for stable fixed points 
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labeling rows (columns) of these matrices as +, — from left to right (from up to down), 
by (|TBp, QT7p and QTB|) we have 

M ( a 3 + a 2 x~ / a 2 + \ _ . . 

y x y ~ 1 (a 2 + aix 1 ) x" / ^ 1 (ai + a x 7 ) J ' 

A = ( i ^ ) and A ' = ( -I ^ ) ■ < 20 > 



Using (|lq) - (p0|) we can rewrite (]15[) as follows 

7Y(.4'M") 

" Tr(AM n ) • (21) 

Using ( |i~3|) it is easy to show that matrices M (|T9|) and A fl20|) commute. Thus, they can 
be diagonalized simultaneously. After performing this diagonalization we have 

(SA)=(^) 2 + 77 ^(mY, (22) 



i + xij {i + xiy \Xi(x) 

where 

Xi(x) = a 3 + 2a 2 x y + aix 27 , (23) 
A 2 (x) = x 7 " 1 (clqx 1 - a 2 x + a x (l-x 7+1 )) (24) 

are the eigenvalues of matrix M. Note that, for stable fixed points (Si) = (S n ) = m. Thus, 
using formula flllj) for m we can reveal the following formula for correlation functions 

, , Ax q f\2(x)Y 

9(n) = M (w) (25) 

where Xi(x) and A 2 (x) are defined in (55), (E3). 

One can prove easily that the correlation function decreases exponentially with the 
distance n independently of x, which is a consequence of the infinite dimension of the 
Husimi lattice. 

For a period doubling bifurcation of the mapping (|TT|) the matrices corresponding to 
the M and A matrices and depending on the x\ and x 2 points of the two cycle {xi,x 2 } 
do not commute, so the procedure applied above is unusable for this case. 



4 Connection with the two- vertex and nn interaction 
Ising models on the Bethe lattice 

In the previous section we have derived a closed formula for correlation functions on the 
Husimi lattice. Now we shall show that for q = 2 our model coincides with the two- vertex 
model on the Bethe lattice with q = 3 . First we recall the definition of the model. Let us 
consider the Bethe lattice with iV shells and the coordination number q = 3. The central 
point O and the adjacent edges are considered as the Oth shell. Each edge, including 
the free ends on the surface of the tree, can be in one of the two distinct spin states 

5 G {+, — }. To each vertex one ascribe the statistical weight it>(Si, S 2 , S3) depending on 
the number of incident edges in the (+)-state. Let denote the corresponding statistical 
weight of the vertex with i incident edges in the state +. Using notation (|3|) and Fig. 2 
it is easy to see now that our models coincide. 
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Fig. 2 The Husimi and Bethe lattices with q = 2 and q = 3 respectively. 



In the following we shall show that our model with only nn interaction ( J 3 = h = 
0) and q = 3 transforms to the nn interaction Ising model on the Bethe lattice with 
q = 3. In this case we shall use widely known star-triangle transformation (or the Yang 
-Baxter equations) [|TJ to transfer our model to the Bethe lattice one. The star-triangle 
transformation replaces a star consisting of three spins interacting with a central spin, by 
a triangle of spins interacting with one another and vice versa (Fig.3). 



Si 



L 




S3 S2 S3 S2 

Fig.3 The star-triangle transformation. 



Let us consider the simplest case in which the interactions of all directions are equal. 
Then summing over the central spin in the partition function of a star and equaling it to 
the partition function of a triangle with some coefficient R, we have 

ex P L (S0S1 + S S 2 + SqS 3 ) = R exp J 2 (S1S2 + S1S3 + S 2 S 3 ) 

So 

where L and J 2 are the pair interaction parameters of the Bethe and Husimi lattices 
respectively. After performing the summation over So and substituting Si = ±1, i — 1,2,3 
we obtain the Yang-Baxter equations in the form 

2cosh(3L) = i?exp(3J 2 ), (26) 

2cosh(L) = J Rexp(-J 2 ), (27) 

From notation @ it's follows 

do = (13 = exp(3J 2 ) and at = a 2 = exp(— J 2 ). (28) 
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Due to summation over spins located inside of triangles in star-triangular transformation 
one can obtain 

9H(n) = #s(2n), 

where gH{n) and gB{2n) are the correlation functions of the Husimi and Bethe lattices 
respectively. Thus, for correlation functions of the nn interaction Ising model on the Bethe 
lattice we have 

4a; 3 
where 

, . a;(exp(2L) — x) 
r(X) = x 2 + exp(2L) ' 

Here we used the fact that x is a stable fixed point of the corresponding recurrence 
relation on the Bethe lattice and it coincides with Husimi's lattice one. Here also one can 
prove that, independently of x, the correlation function decreases exponentially with the 
distance n, which is a consequence of the infinite dimension of the Bethe lattice. 



9B(n) = ; ,, 2 r n {x) 



From (p6|) , ( f27|) after simple algebra one can obtain 

sinh2Lsinh2J 2 = k^ 1 (29) 

where 



sinh d 2L 



2 (cosh 3L cosh 3 L 



(30) 



First the star-triangle transformation was used for triangle and honeycomb lattices [T|]. It 
was obtained that for triangle and honeycomb lattices the condition of the critical point 
is k = 1. For the ferromagnetic Ising model on the Bethe lattice it is well known that the 
critical point is 

L r = - In 



2 q-2 

where q is the coordination number of the lattice. From ( |30"D follows that for the case 
q = 3 parameter k is unequal to 1 (k — g^)- From ( p9"D one can obtain the critical point 
of the ferromagnetic Ising model on the Husimi lattice ( J 2 )c = \ In f (the last result can 
be derived in a straightforward manner, see Sect. 2). 



5 Concluding remarks 

In this paper we have derived a closed formula for correlation functions of the multisite 
interaction Ising model in the presence of an external magnetic field. We showed that the 
correlation function decreases exponentially with the distance which indicates no existence 
of long range correlation at any finite temperature this is a consequence of the infinite 
dimension of the Husimi lattice in Euclidean space. This result is in good agreement with 
the results obtained in papers || ||. All our obtained formulae contain only statistical 
weights dj, which depend on the spins located in the vertexes of triangles. Hence, our 
formulae can be used also for other types of Hamiltonians where nn (nearest neighbour) 
and three-site interactions are present, for example Potts model, etc. 
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The star-triangle relation (the Yang-Baxter equations) maps a high-temperature model 
on the Bethe lattice with q = 3 to the high-temperature model on the Husimi cactus. 
Together with nonanalytic behavior of the recurrence relations of the Bethe and Husimi 
lattices one can use them to obtain the critical exponents and correlation functions too. In 
conclusion, let us note that the results deduced on the Husimi and Bethe-like lattices are in 
good agreement with results obtained on planar lattices, and according to Gujrati, Bethe 
or Bethe-like (Husimi cactus) calculations are more reliable than conventional mean-field 
calculations |22| . 
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